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I. INTRODUCTION 



The mean-field approach to the phenomenon of superconductivity introduced by 
Bardeen, Cooper, and Schrieffer (bcs) allows for a simple and elegant treatment of 
pairing correlations in nuclei The main drawback of the BCS method is that its wave 

function is not an eigenstate of the particle number operator. The accuracy of the BCS 
approximation is satisfactory if the pairing interaction strength is strong or the number of 
particles is very large These conditions are not satisfied in nuclei. Indeed, the critical 
value of the effective pairing strength, Gcrit, above which the static gap exists, is inversely 
proportional to the single-particle level density around the Fermi level [Q, i.e., it becomes 
very large around subshell closures. 

It was demonstrated by Lipkin |^ that the effect of the nucleon number fluctuation can 
be suppressed by using a model Hamiltonian H — XiN — \2N'^ instead of the conventional 
Routhian H — XN, where H is the original Hamiltonian (involving pairing interaction) and 
N is the nucleon number operator. The approach by Lipkin was then developed by Nogami 
and his collaborators []7|-|ri|l. The important feature of the Lipkin-Nogami (ln) method is 
that (i) there always exists a nontrivial (superfluid) solution regardless of the strength of 
pairing force and (ii) the LN wave function has a similar form to that of the BCS method, 
thus allowing a simple interpretation of excited states in terms of quasiparticles. 

In this respect it seems to be of a considerable interest to extend the LN method to 
the case of rotation, where the short-ranged attractive pairing force plays a significant role. 
For instance, at low spins the pairing correlations tend to significantly reduce the nuclear 
moment of inertia as compared to the rigid-body estimate. On the other hand, it is well 
known that at very high spins many nuclei behave as macroscopic rotors, i.e. their moments 
of inertia are fairly constant and close to their rigid body values. More examples illustrating 



the importance of pairing correlations at high spins can be found in the review in Ref. |T2 

In this connection further studies of the Mottelson-Valatin effect [|l3l (phase transition 
from superfluid to normal state in rotating nuclei) are of great importance [p!^p!2| . Since the 
BCS method provides a rather poor description of the pairing phase transition region, one 
expects that the LN method would be a powerful tool allowing for a better description of 
pairing correlations without losing the simplicity of the rotating independent quasiparticle 
picture. 

In Sec. 2 the cranked LN (lnc) equations are derived. The equation for the particle 
number fluctuation is explicitly written in terms of the single-particle and pair densities. 
This will be useful when adopting the LN approach to the general Hartree-Fock-Bogolyubov 
(hfb) treatment. The method is applied in Sec. 3 to the exactly solvable two-level cranking 
model with pairing. The exact solutions are compared to those obtained within the standard 
mean field approach, particle number projection method, LNC treatment, and its particle- 
number projected version. The conclusions are presented in Sec. 4. 



II. THE METHOD 
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A. Cranked Lipkin-Nogami equations 



Let us consider the BCS Hamiltonian which contains a single-particle Hamiltonian, H^p, 
and a seniority pairing (monopole, state independent) interaction: 

H = Hsp + -^pair = ^kCitak " G ^ ala^ttjai, (2.1) 

k k,l>0 

where is the single-particle energy, G is the pairing strength, and \k) = T\k). 

In order to investigate the pairing interaction in rotating system we consider the cranking 
Hamiltonian (Routhian): 

H"^ = Hsp + i/pair — ^^jx, (2.2) 

where denotes the component of the total nucleonic angular momentum on the rotational 
axis (here: x-axis). This axis is assumed to be fixed in space (only one- dimensional rotation 
is considered) and the angular velocity of rotation, uj, is supposed to be constant. 

In the standard rotating BCS (rbcs) approach the expectation value of the Routhian 
( p.2|) , E'^ = {H'^), is minimized in the product state of independent quasiparticles defined 
through the Bogolyubov transformation: 

' (2 3) 

I 

where af{ai) is the quasiparticle creation (annihilation) operator. The RBCS wave function 
represents a mixture of states with different numbers of particles. Consequently, in order to 
account for the fluctuations, the particle number should be projected before variation. This 
can be done in a good approximation by means of the LN method outlined below. 
Let us assume that the state \ipn) is the quasiparticle vacuum, i.e., 

= 0. (2.4) 

The index n stands for the average number of particles in the state {ipn), determined by 
means of the particle number equation 

{^n\N\iJn) = n, (2.5) 
where N = ^a^a^ is the particle number operator. The state \ipn) can be expanded in 

k 

eigenstates of the particle number operator, 

l^n) = ^C„,„j0„o), (2.6) 
no 

iV|0„„)=no|0„o). (2.7) 

The total Routhian ( p.2|) commutes with the particle number operator, [H'^jN] = 0. By 
expanding E'^ in terms of n 



3 



{(t)n\H''\(t)n) = AH = Ao + Xin + X^n^ + ... (2.8) 
the following relation is obtained: 



I^"l0no) = (^n|i^"|^n) " Ai ( | iV| V'„) - U^) 

- A2((V^„|iVVn)-n^)-... (2.9) 

The important feature of the above expression is that it yields the expectation value of 
in the projected RBCS state, |0no); in terms of the expectation value of Routhian in the RBCS 
state, \ipn)- This implies that, knowing coefficients Aj, one can minimize the right hand side 
of Eq. ( p.9|) instead of minimizing explicitly the expectation value of in the projected 
RBCS state. 

So far the considerations are exact. In the next step, however, the expansion (|2.9|) is 
truncated by retaining the first m terms; the coefficients Ai,A2,...,Am are then calculated 
from the following set of linear equations: 

{^n\K^N\^^) = {i,n\K^\^n){'4^n\N\'ipn) 

{ijn\k-N^\i^n) = mk^\i^n)mN^\i'n) (3 10) 

where = — XiN — A2iV^ — ... — XmN"^. The m = 1 case, together with the constraint 
( ^TSD , is equivalent to the standard RBCS method. The case of m=2 is discussed below. The 
generalization to higher orders is straightforward although algebraic manipulations become 
tedious. Since the right hand side of equation (|2.8| ) has to be finite, the i-th order correction 
to Routhian, which is proportional to Aj is of the order of n~' . Therefore, when one deals 
with a system with large number of particles it seems reasonable to neglect the third order 
term (the mass formula derived from the seniority scheme is actually quadratic in n). 

In the LN method the expectation value of K''^ is minimized assuming that coefficients 
Ai and A2 are constant (the simultaneous variation of A's would lead to more complicated 
equations 0): 

= 0. (2.11) 

The Lagrange multipliers A and u are determined by fixing expectation values of the particle 
number and angular momentum, respectively, 

(tPnim^n) = n, (2.12) 

(^nlilV^n) =/. (2.13) 

It is convenient to express the results through the density matrices p and u and the pairing 
tensor x, which are defined by means of the transformation matrices A and B: 

Pki = J^iBl^Bii, 

Uki = Ei Al,Aii = 5h - pkh (2.14) 
Xki = Ei^K-Bfej. 
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Expressed in terms of quasiparticles the operator takes the form 



~ -"-00 ^ -"-11 -"-20 ^ -"-22 -"-31 ^ 



40; 



where 



-^00 = XI [^fc^fc' ~ pki 



k,l 



A2 



~ oG* H PklPM 



k,l 



\in - Aan^ - 2A2 Pfc/W^, 



(2.15) 



(2.16) 



k,l 



+ X2Sij \afaj, 



(2.17) 



i,j (. k,l 



(2.18) 



The terms K22, K^i, K'^o represent the residual interaction between quasiparticles and are 
neglected in this approximation. In the above relations the following quantities are intro- 
duced: 



^ki = ekSki - Gsign{k)sign{l)p],i + AX2PI1, 



(2.19) 



A = Ai + 2A2(n+ 1), 



Am = -5fe,sign(A;)A = -G6kisign{k) ^ Xkk- 

k>0 



Condition (p.ll|) leads to the hfb equations: 

I 



where 



^kl = ^kl - UjUx) 



St = Et - \, 



kh 



(2.20) 
(2.21) 



(2.22) 



(2.23) 
(2.24) 
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The difference between the usual hfb equations and the above ones consists in the appear- 
ance of the parameter A2 which should be determined selfconsistently from Eq. (|2.10|) . It is 
important to note, that the eigenvalues of the ln+hfb equations (|2.22|), , are related to 



quasiparticle Routhians E"^ through relation ( p.24|) . Consequently, special care should be 
taken when interpreting the results using the standard technique of quasiparticle diagrams 
of the cranked shell model. 

The total LNC energy of the quasiparticle vacuum is given by 



E- 



LNC 



1 

X! ^kPkk 77 ~ X! PklPkl 

k ^ k,l 

'^^2'^ PklUkh 
k,l 



(2.25) 



where the term proportional to A2 represents the nucleon number fluctuation correction. 

The presence of selfconsistent symmetries very often facilitates the calculations. One 
such symmetry, important in the context of cranking model and high spins, is the signature 
symmetry, i.e., the symmetry with respect to the rotation of the system by 180° around the 
X-axis. The single-particle states with good signature r are related to the original fermionic 
basis by the so-called Goodman transformation Wl 



\K,r 



1 

1 
72 



-\k) + (-1)"^-2|A;) 

+ \k) 



(2.26) 



where is the projection of the single-particle angular momentum on the symmetry axis. 
Application of transformation ( p.26| ) leads to immediate selection rules for the coefficients 
of the Bogolyubov transformation. 



A 



KL 



B 



KL 



0, 



(2.27) 



and, consequently, for the matrix elements of the ln-|-hfb Routhian [15|. In the following. 



it will be assumed that the system is invariant with respect to the signature symmetry. 



B. Calculation of Xo 



In order to compute A2 one can use the set of equations (|2.1CI| ) which in the case of LNC 
is reduced to the two relations: 



(2.28) 



Because of the requirement = 0, the first equation is satisfied automatically and the 
second one is simplified to 



(0|i^'^|4)(4|iV2|0) = 0, 



(2.29) 
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where |0) denotes the quasiparticle vacuum and |4)(4| is the projection operator on the 

four-quasiparticle space. The matrix elements of and A^^ that appear in ( |2.29| ) are given 
by 

^40 = H I T ^KP^ B^p^ Aip^ Alp^ 

K,L,P I ^ 

+ ^2BKPiAKP2BLPs^LPi papiap2ttP3ttP4 



+ h.c, (2.30) 



^40 - 51 ^ArQi-^JVQ2^MQ3-^M(34"Qi"Q2"Q3"Q4 
M,N,Q 

+ h.c., (2.31) 



where we have aheady apphed the Goodman transformation to states with well defined 
signature. By means of relations ( p.3CI| ), ( p.31| ), and ( p.29| ) one can now obtain the expression 
for A2: 

„ {XKdP*KL + PKl)} {XKli'^kL+UKL)} -'^Y.^XX^)LK{XX^)Lk 

X, = ^^^ ^i^^^^^ ^ . (2.32) 

4 [Tr(xx+)] -2Tr(xx+XX+) 

It is easy to show that for 00 = the above relation reduces to the well known result (see 
Ref. [|TT|): 



X _ G_k>0 fc>0 fc>0 /r, oo\ 

fc>0 fc>0 



Relation ( ^.32| ) together with equations ( |2.22| ) completes the derivation of the LNC equations. 



III. RESULTS FOR THE R(5) MODEL 

In order to examine the accuracy of the LNC approximation we consider the two-level 
Krumlinde-Szymahski R(5) model (see |T6|-p!^, p!^ , p!9[] ) . The Hilbert space of this model 



consists of Q j=| multiplets. As shown in Fig. |T], the single-particle levels are split by 
the deformation of the average nuclear potential. The single-particle splitting is 2e, i.e., 
the energy of the upper levels (labelled as |l) = |m = 3/2) and |1)=T|1)) is +e and that of 
the lower levels (labelled as |2) = |m = 1/2) and |2)=T|2)) is — e. We begin our analysis by 
considering the half-filled (symmetric) system, i.e., the number of particles is equal to n=2Q. 
In the particular version of the two-level model considered in this paper, the Coriolis coupling 
between the lower levels is neglected, i.e., {jx)22 = 0. The two-body pairing interaction is 
assumed to be of the monopole type. This simple model contains the essential physical 
features of the nuclear structure relevant to the interplay between pairing and rotational 
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motion. Its main advantage is that it can be solved exactly using the Lie algebra associated 
with the symmetry group -R(5) [|17 . 

In the standard notation of representing HFB equations in doubled dimensions 0, the 
single-particle field v has the form of a 4 x 4 matrix: 



V 



K 





(3.1) 



where 



-ijj 



Gpii + 4A2Pn, - Gpi2 + 4A2P12 
Gp-2i + 4A2P21 . e - Gp-22 + 4A2p;2 



(3.2) 



-e - Gpii + 4A2Pn, ^ - <^Pi2 

UJ - Gp2l + 4A2P2i, ^ ~ ^P22 



■ 4A2P^2 

4A2Pi2 



In the above relation it was assumed that fl4,l£ 



{3x)2l = {jx)l2 = -{jx)21 

Similarly, the pairing field is given by 



A = -Gnixn + X22) 



12 



1. 



(3.3) 



(3.4) 



" 





1 


0" 











1 


-1 














-1 









(3.5) 



The LNC equations are solved iteratively with respect to the density matrix and pairing ten- 
sor. The initial values can be taken from the standard RBCS equations. Then the coefficient 
A2 is found and the matrices A and B are computed from Eq. (|2.22| ). This gives a new 
approximation for the density matrix and pairing tensor. One should choose properly A (the 
chemical potential) at every iteration step to satisfy the relation n=Tr(p). The above pro- 
cedure is continued until the convergence is achieved. For more complicated Hamiltonians 
it is suggested to use the so-called gradient method (see Refs. [pOH 2^). 

Having found the Bogolyubov transformation matrices A and B, one can calculate p and 
M, Eq. ( p.l4|) , and then the total energy, 



tot 



ef2(pii + pii - P22 - P22) 
1 

77 - i^GVt {pMpki) - 2A2l^ (pklUkl), 

^ ^ k,l>0 k,l>0 



(3.6) 



the pairing energy. 



pair 



= -£^tot 

+ 2ne 



E, 



unpair 



E, 



tot 



2?^ 



2n 



(3.7) 
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the pairing potential (energy gap), 



A = + X22] 



(3.8) 



the total angular momentum, 



I = Tr(j^p) = 2n{pi2 - P12). 



(3.9) 



In the light of the recent results [^, it is advantageous to carry out the particle num- 
ber projection after Lipkin-Nogami (lnc+pnp). This can easily be done by following the 
formalism of Ref. fl^. The particle- number projected energy and angular momentum are 
given by 



where 



Ei! = -2Qe 



P12 Pn-iiO 
P{+) - P(-) PniO 



^sp 



is the single-particle energy, and 

2 k'^p{+)p{-) /2p(+)p(_) P/j_i(0 



(3.10) 



(3.11) 



(3.12) 



4\ = -2Gn' 
and 



+ [p(+) - p(_)] 



M+)PI) Pn-iiO 



(3.13) 



1 - 



(l]-l)[p(+)-p(_)P PniO 



+ 4:P 



12 



\ _ , Pn-i(0 1 1 

(3.14) 



are two contributions to the pairing energy (cf. Eqs. (A. 7) and (A. 8) of Ref. ^^), and 

k = p22-pu, (3.15) 



Pi±) 



(l±V'F + 4py, 



PW + pI) 

Pi+) - P{~) ' 



(3.16) 



(3.17) 



In the above equations Pn{x) is the Legendre polynomial of the n-th order. It is easy to see, 
that in the limit of very weak pairing (^ ^1, P„(l) = l) the total pairing energy given by 
Eqs. (|3.13|) and (|3.14|) becomes equal to -Eunpair- On the other hand, if pairing is very strong 
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^oo, Pn-i/Pn n/{2n - 1)^ ^ P^_^/Pn -> n{n - l)/(2n - 1)^ \ the pairing energy 
approaches the hmit of the seniority model, — + 1). 

We have performed calculations within the symmetric variant of the R(5) model for e=l, 
f2=20 (i.e., the half-filled symmetric system with n=40 particles), and for three values of 
the pairing strength, (7=0.015, 6^=0.065, and (7=0.1 ||2^. Without rotation, the mean- 
field (bcs) solution undergoes a transition to the paired regime at the critical strength 
(?=2e/(2r2 — 1)~0.051. Therefore, the intermediate value of (^=0.065 corresponds to the 
phase transition region, and allows us to study the destructive role of rotation on pairing 
correlations, while the other two values of G represent the weak and strong pairing limits, 
respectively. 

The results of the BCS method are in Figs. |^-^ denoted by the full triangles and full 
circles for RBCS and RFBCS, respectively. The latter results are obtained by the variation 
after projecting the good particle number component of the RBCS state. The results based 
on the LNC method are denoted by open symbols. The open triangles and open circles 
represent the LNC and lnc+pnp (exact particle number projection of LNC solutions) results, 
respectively. The exact results are denoted by the water-wheel symbols. 

Calculations have been performed for lo ranging from to 1.2 [^]. For each cu, the energy 
and angular momentum have been determined from Eqs. ( |3.6D and ( |3.9D , or Eqs. ( |3.10| ) and 
( p. 11 ), and the pairing energies have been computed by means of Eq. ( p.7| ). In this way, the 
plots of pairing energy versus angular momentum have been constructed and are shown in 
Fig. |. 

For the weak pairing strength. Fig. |^a, the RBCS method gives only the unpaired solution 
for all spin values [0. Although the solutions of the LNC method contain some pairing 
correlations, the approximate formula for the energy with the corrective A2 term, Eq. ( p.6|) , 
gives the pairing energy much too small. On the other hand, when the LNC solutions are 
projected on the good particle number one obtains a fair qualitative agreement with the 
exact results (the maximum relative difference is of the order of 20%). This illustrates the 
fact that the LNC wave function describes resonable well the pairing correlations in a weak 
pairing limit even if the average value of the auxiliary Hamiltonian is not a very good 
approximation to the exact energy of the system. 

As seen in Fig. |^, the RFBCS results provide an excellent approximation to the exact val- 
ues. This indicates, that there is still some room for improvements of the LNC+PNP method. 
From Ref. p3| it follows that the difficulties of the LN method in describing the half-filled 
(n=2f2) system in the weak pairing limit stem from the fact that its exact ground-state 
energy cannot be approximated by a second-order expansion in n centered at no=2fi. How- 
ever, the parabolic expansion works very well for the ground-state energies of asymmetric 
systems with n7^2f2. Based on this result, a useful two-step procedure is suggested. Firstly, 
the LN or LNC equations are solved for the system with no=2f2 ± 2. In the second step, 
the ?7,=2n component is projected out from the resulting LNC wave function. Such a hybrid 
method relies on extrapolating the n=2Q solution from either those for n <2Q or those for 
n >2Q, instead of interpolating between solutions for n <2Q and those for n >2Q. Results 
of LN+PNP calculations for no=2fi+2=42 are presented in Fig. § and agree remarkably well 
with those using the RFBCS method and with the exact values. 

The weak-pairing regime is realised in nuclei around shell or subshell closures where 
e^Gfl. Therefore, the hybrid method described above can be useful for studying properties 
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of, e.g., spherical magic nuclei at low spin or superdeformed magic nuclei (^^^Dy or ^^^Hg) 
at high spin. 

For the intermediate pairing strength. Fig. ^b, the RBCS method yields the unpaired 
solution above a certain critical angular momentum. Such a sharp transition is not present 
in the exact results. Although the low-spin LNC results are in much better agreement with 
the exact solutions than those of the RBCS method, this method yields too small pairing 
energy at higher spins. Not surprizing, the pairing energies of the pnp+lnc method agree 
quite well with those obtained by means of the RFBCS treatment. The latter ones describes 
fairly well the exact values, while the remaining small difference, which cannot be accounted 
for by using the BCS-type wave function, illustrates the presence of large correlations (due 
to the quasiparticle interaction) at the phase transition region. 

For the strong pairing limit. Fig. the RBCS method is the only one which fails to 
describe exact results. All other methods give very good agreement at low and high angular 
momenta and leave a gap of missing solutions around J=18, i.e., in the region where the 
adiabatic cranking approximation breaks down (the solutions in the phase transition region 
correspond to the maximum of the total Routhian rather than to the minimum l]!^). Inter- 
estingly, the region of instability is very large for the RFBCS method, where it extends down 
to 1=10, while the LNC method is able to follow solutions up to J=16. The particle number 
projection does not lead to a significant improvement here. 

Fig. ^ shows the pairing delta Abcs as a function of uj. This quantity characterizes 
the BCS states used in RBCS, LNC, and RFBCS methods and defines the BCS occupation 
probabilities (for the LN occupation probabilities, see, e.g., Ref. [^). In fact, in the methods 
employing the concept of pnp, Abcs is not related to any particular observable; it should not 
be understood as the energy gap, but rather as a variational parameter (see Ref. fl^). Fig. 



^ shows Abcs together with Aexact='\/—G'( impair) where the average value is calculated with 
respect to the exact ground-state wave function. By definition, Agxact is a direct measure of 
pairing correlations and it should reflect Abcs in the limit of large pairing. Note that the 
Abcs parameter for lnc-I-pnp is the same as for LNC, and, therefore, is not shown in Fig. 
^. It is seen in Fig. ||a that in the weak pairing limit the two quantities Agxact and Abcs 
differ very much even if the RFBCS pairing energy is rather close to exact values. Fig. ^. In 
the intermediate and strong pairing limits the Abcs parameter of RBCS goes to zero at the 
critical angular momentum, which illustrates the Mottelson-Valatin phase transition. The 
exact results do not show such a sharp transition, and the Abcs parameters of the LNC and 
RFBCS method qualitatively reproduce the exact values. 

The angular momentum alignment is illustrated in Fig. | where the kinematical moment 
of inertia, 

J(') = — , (3.18) 

is shown as a function of rotational frequency uji. At this point it should be stressed that 
the rotational frequency, cj/, obtained from the canonical relation 

dE 

= -Jf (3.19) 

is equivalent to the cranking-model frequency uj only for the exact eigenstates of a model 
cranking Hamiltonian. Consequently, the relation 
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(3.20) 



holds exactly for the exact solutions of the R(5) model, and also for the RBCS variant (the 
solution is an eigenstate of the independent-quasiparticle Hamiltonian) . On the other hand, 
the relation ( |3.20D does not hold for the approaches based on the pnp treatment, since the 
resulting states are determined from the restricted variational principle. Consequently, the 
rotational frequency in Fig. ^ is determined using Eq. ( 3.19|) . The moment of inertia, Eq. 
( p.l8| ), illustrates the Mottelson-Valatin phase transition. It can be seen that for all pairing 
strengths the LNC method correctly describes this transition, both with and without the 
subsequent particle number projection. 

We have also performed calculations for the asymmetric system, n^2VL. Here, the static 
RBCS solution is always present, independently of the value of the pairing strength and 
the results and conclusions are very similar to those for the symmetric system, n=2f2, in the 
strong pairing limit. The results for the n=48 system are shown in Fig. ^ It is seen that 
the LNC method provides an excellent approximation to the pairing energy, even without a 
subsequent particle number projection. 



IV. CONCLUSIONS 

We have studied the pairing correlations in rotating nuclei using the cranked Lipkin- 
Nogami LNC method which is based on employing the auxiliary Hamiltonian = H — 
XiN — X2N^ — ujjx, where the parameters Ai and A2 are chosen so that the influence of the 
nucleon number fluctuation is strongly reduced. One should emphasize the simplicity of the 
LNC approximation, especially when compared with more sophisticated projection methods. 
In practice, the LNC method is a simple extension of the usual RBCS treatment. 

Good accuracy was obtained for the ground state energy, particularly in the case of strong 
pairing interaction. It means that the method suppresses correctly the "dangerous" (particle- 
number violating) part of the quasiparticle interaction. The weakness of the LNC method 
in the weak pairing limit can be overcome by performing the projection after variation. 
Therefore, the lnc+pnp approach can provide us with a fair description of rotating nuclei 
near shell and subshell closures. 

Another welcome feature of the LNC method is that it provides us with a very good 
description of the pairing phase transition region, regardless of the strength of pairing inter- 
action. Note that using the LNC method one can obtain non-trivial solutions even for very 
fast rotation where the RBCS method breaks down. 
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FIGURES 



FIG. 1. Energy levels for the R(5) model. The above pattern is repeated times. 

FIG. 2. Pairing energy versus angular momentum for three values of pairing strength, G=0.015 
(a), G=0.065 (b), and G=0.1 (c). 

FIG. 3. Pairing energy versus angular momentum for the weak pairing strength, G=0.015. 
The LNC+PNP results calculated for n=2f2=40 from the no=42 solutions are also shown. 

FIG. 4. Pairing delta A versus rotational frequency ui for three values of pairing strength, 
G=0.015 (a), G=0.065 (b), and G=0.1 (c). 

FIG. 5. Kinematical moment of inertia versus rotational frequency for three values of 
pairing strength, G=0.015 (a), G=0.065 (b), and G=0.1 (c). 

FIG. 6. Pairing energy versus angular momentum for the asymmetric variant of the R(5) 
model with n=48 and for three values of pairing strength, G=0.015 (a), G=0.065 (b), and G=0.1 
(c). Only the exact, RBCS and LNC results are displayed. 
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